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We examine the Q 2 evolution of gluon polarization in polarized nucleons. As is well known, the 
evolution of a s AG(Q 2 ) is negligible for typical momentum transfer variations found in experimental 
deep inelastic scattering. As a s increases, however, the leading nonzero term in the evolution 
equation for the singlet first moment reduces the magnitude of the gluon spin. At low Q 2 the term 
a s AG can vanish, and ultimately become negative. Thus, low energy model calculations yielding 
negative AG are not necessarily in conflict with experimental evidence for positive gluon polarization 
at high Q 2 . 
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I. INTRODUCTION 



Polarized deep inelastic scattering (DIS) has proven to be an important new tool for studying the structure of the 
nucleoli. Through sum rules, the experiments have provided us with values of various operator matrix elements; for 
instance, the singlet sum rule indicates that quark spin accounts for only a small fraction (s=s 1/3) of the nucleon's 
total spin p]|^]. One possible resolution of this apparently surprising result is that the gluon fields in the nucleon 
bound state may carry a significant fraction of the nucleon's spin. Orbital angular momentum of quarks and gluons 
may also be significant. 

Recently, various authors |3|JJ] have attempted to make estimates of the size of the gluon contribution to the spin 
of the nucleon, 



12 



AG{Q 2 ) = dxAg{x,Q 2 ) = dx(g^x, Q 2 ) - 9l (x, Q 2 )). (1) 
Jo Jo 

In at least one case ||, this quantity was evaluated in low energy models, and turned out to be negative. This is 
in apparent conflict with extractions, primarily at DIS scales, which find AG positive (^,^|. Indeed, to leading order 
(LO) in perturbative QCD, a s {Q 2 )AG{Q 2 ) is renormalization group invariant, thus the claim has been made that at 
least the sign of AG should be a reliable prediction of the models 

It is well known j^] that the singlet part of the first moment of the spin structure function g\ has an anomalous Q 2 
evolution. It is also well known that the leading term in a s in the axial anomalous dimension vanishes, and it 
is for this reason that many authors dismiss this Q 2 evolution as insignificant. Roughly speaking, corrections to the 
singlet first moment arising from the anomalous dimension can be argued to behave like a s log Q 2 , and hence appear 
approximately Q 2 independent. This is, however, not precisely true at higher order in a s . In an earlier paper, j^] 
it was argued that comparing momentum sum rules from unpolarized electroproduction and the spin sum rule for 
<7i, including their Q 2 evolution, showed that DIS spin measurements are consistent with a low energy valence quark 
picture, where the valence quarks carry a substantial part of the spin of the proton. In a later paper, |l0[ , those 
calculations were extended to consistently include the next higher order QCD corrections in leading twist. Here, we 
will explicitly examine the effects on gluon spin. We evaluate the leading non-zero QCD corrections to the gluon spin 
evolution, and find that the sign of the first correction term forces one in the direction of smaller (more negative) 
a s AG at lower Q 2 . Application of these perturbative formulae show that a sign change in a s AG thus need not be 
such a surprise. 

There are several caveats to this analysis which must be made at the outset. First, we apply perturbative formulae 
down to a regime where they cannot rigorously apply (a s /iT « 1) and thus do not expect to have quantitative 
predictive power. However, the sign change in AG is significant, and if the perturbative result goes in this direction, 
there is certainly no reason why the full nonperturbative result might not behave similarly. Furthermore, one can 



point to the poorly understood successes of the one-gluon exchange mechanism in the hadron mass spectrum 1 1 1 , or 
of the evolution of unpolarized structure functions j^] , which provide some faith in at least the qualitative predictions 
of this method. 

The second issue is one of scheme dependence. The decomposition of the singlet sum into terms involving AS and 
AG(Q 2 ) depends on renormalization and factorization scheme choice, and is gauge dependent Indeed, there is no 
unique way to define quark and gluon distributions beyond LO. When in a perturbative QCD regime, we work in the 



Adlcr-Bardeen (AB) factorization scheme |3|,|3 10 13J, with MS renormalization conventions, which defines AG and 
its contribution to the singlet first moment as well as its Q 2 evolution to LO. However, if one wants to work beyond 
leading order, in particular to evolve to low scales to compare with model predictions, one must know in addition the 
connection with the choice of scheme (and thus the underlying operator matrix element definition) of e.g. AG in the 
low energy model calculations. Thus, we really cannot make any direct quantitative comparisons with the numerical 
values extracted in quark model calculations. We can, nevertheless, unambiguously examine the leading evolution of 
the gluon spin downwards in Q 2 given our choice of scheme, which shows a clear trend towards zero and eventually 
negative values. Ultimately, one must understand both the renormalization scale and scheme dependence of the quark 
model calculations in order to decide if the contribution of gluon spin is quantitatively compatible at low and high 
energy scales. 

II. FORMALISM 

Working at leading twist, the contributions to the hadronic tensors involved in polarized DIS, via the operator 
product expansion, are due to singlet and non-singlet axial currents. The first moment of the polarized spin structure 
function gi(x, Q 2 ) is given by 
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if n > = fdxg^ = \c s (Q 2 )a Q {Q 2 ) + \c NS (Q 2 )(±a 3 /2 - 



as/6) (2) 



where C S (NS) are (non) singlet coefficient functions calculable in perturbative QCD which have been worked out to 
second (third) order in a s Both are given by 1 — a s {Q 2 )/i^ to leading order. 03 and as are renormalization 

group invariants, and can be obtained (modulo issues of flavor SU(3) symmetry breaking fli]) from low energy /?-decay 
constants. 

The singlet matrix element is obtained from 

n f 

ao(Q 2 )s^ = (p,s\J^\p,s) = (p,s\^2qa^ 5 qi\p,s), (3) 

i=i 

which in a naive quark model picture is just the total spin of the quarks. Experimental results are generally given 
for values of ao(Q 2 ) with the coefficient function C S (Q 2 ) factored out, as shown above, but because of the UQ) 
anomaly in QCD, ao is itself not scale independent. The QCD evolution of ao has been worked out to NNLO [p[, 
and is given below. For our purposes, we also wish to consider a decomposition of ao into separated quark and 
gluon contributions. This is a renormalization scheme-dependent, and gauge dependent, separation. Some authors 
e.g. choose a so called gauge- invariant scheme in which the gluons by definition do not contribute to ao E3]. In 
the Adler-Bardeen factorization scheme [^|Jl3|, sometimes called a chirally invariant scheme, the anomalous gluon 
contribution is explicitly separated out as follows: 

a (Q 2 ) = AH- fa s (Q 2 )AG(Q 2 )/2ir (4) 

where f is the number of active flavors at this Q 2 . The evolution of the full singlet first moment, ao, is known to next 
to leading order (NLO) in a s 

a (Q 2 ) = a {Ql) exp - / da'^-{ (5) 

where j s (a s ) = 7^ (as/An) 2 + 7^ (a s /4n) 3 + ■ ■ ■ is the singlet anomalous dimension, with 7^ = 16/, and 7^ = 
(944//3 — 32/ 2 /9) with our choice of renormalization conventions, and (3 is the beta function of QCD, 

2 3 

(3{a s ) = da s /dr = -fcSL-fr-2!— + ,.. (6) 

47T l07T z 

where r = ln(Q 2 / ' k 2 QC D ) , and fa = 11- 2//3, ft = 102 - 38//3 U|. To leading non-zero order, this gives 

a (Q 2 ) » a (Q 2 ) (l+g^ («,(Q 2 ) - a.(Q2)) + 0(a 3 s )\ (7) 

Note that the difference a s (Q 2 ) — a s {Q 2 ] ) appearing in the formulae above is itself of order a 2 , thus naturally the 
evolution of ao vanishes at order a s . The solution for ao including NLO corrections in the singlet anomalous dimension 
and the beta function is |p|,^p[ 

ao(Q 2 ) = ao(Q 2 ) (l + ^(^(Q 2 ) - a s (Q 2 )) + (M__^) ( a 2 s (Q 2 ) a 2 s (Q 2 )) 

'■' ,2 -(a s (Q 2 )-a s (Q 2 )) 2 + 0(ai)), (8) 



128tt 2 /3 2 
truncated at order a 3 .. 

It is important to emphasize that the evolution of ao(<5 2 ) is well defined, and is not itself gauge dependent. It is the 
separation of ao into "spin" and "gluon" terms which is scheme and gauge dependent. Using the separation choice of 
Eq. (Eh , the evolution of these quantities are given by @ [l(J 

d_( AE(r) \ 2 / 7SS 7sr W AS(r) \ 

dr \a s {T)AG(r)llTx ) sK ' ^7rs 7rr J Wt)AG(t)/2tt J [ ' 

To leading order, 7rs = 1/27T 2 , 72s = 0, 7sr = 0, 7rr = —f/2ft 2 |§- The evolution manifestly begins only at second 
order in a s . To leading non-zero order, solving Eq. (|^) directly (or alternatively, combining Eqs. |]with 0) we see 
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a s (Q 2 )AG(Q 2 ) = a s (Ql)AG{Ql) - -f- (a s (Q 2 ) - a s (Q 2 )) a (Q 2 ). (10) 

Of course, to first order in a s , the above equation states that q s AG is a renormalization group invariant, the usual 
result. This last equation already shows qualitatively the results we have claimed; namely, if Q 2 is some large scale 
appropriate to DIS, and Qq is some low scale appropriate to a quark model calculation, the leading nonzero correction 
to a s AG is large and positive, and there is no reason the sign need be preserved between low and high Q 2 scales. In 
the Adler-Bardeen scheme, Eq. does not receive higher order corrections j5), and one can combine Eq. (|I|) with 
Eq. (|J) to examine the NLO modifications to a s AG. 



III. RESULTS 

For dcfinitcncss, we use the following experimental values §: a s (5 GeV 2 ) = .287 ± .02, a (5 GeV 2 ) = 0.37 ± .11. 
The evolution equation, (||) evaluated to LO (NLO) yields 

a (5 GeV 2 )/a (10 GeV 2 ) = 1.008(1.010), (11) 

which is a negligible change, well below the present limits of experimental observability. In this sense, the fact that 
to leading order there is no Q 2 evolution is born out. Leaving the region of small a s , however, yields some nontrivial 
evolution. As discussed above, the separation of clq into AE and AG requires additional experimental and theoretical 
assumptions; recent analyses |2|,|| consistent with the scheme conventions of this work finds AG(5 GeV 2 ) ranging 
from 0.8 to 2.6, and thus a s AG(5 GeV 2 ) ranging from 0.2 to 0.75, and AS (which is scale independent here) ranging 
from 0.4 to 0.65, depending on how the analysis and extraction are done. The most naive predictions originally gave 
AE = 1, but many effective low energy quark models are consistent with a lower value, around AE w 0.65 ||. (The 
reduction from unity might come e.g. from the lower components of the relativistic quark spinors, the same source 
which reduces the axial charge in the bag model from 5/3) Using any such prediction in Eq. (|4j), however, assumes a 
partonic interpretation of AE, which is by no means required by QCD, and may be entirely incorrect Jl6[ |. 

Starting from the above experimental values, a s AG will evolve from somewhere in the range (.2 <-> .75) at 5 GeV 2 
down to zero at a scale where a s /ir is in the range (.3 «-> 1.5), i.e. where nonperturbative physics should just begin to 
become significant. Going further down in Q 2 , a s increases and a s AG becomes negative. A smaller assumed value for 
ao (5 GeV 2 ) results in less rapid downward evolution. In Figs. 1-3, we show the result of evolving the singlet moment 
as well as the gluon spin up in a s (down in Q 2 ). If one applies the NLO corrections of Eq. (||) to ao, this makes the 
evolution slightly more rapid. That is, one need not go so far in a s /Tt to get to the "turnover" point (AG = 0, or 
ao ~ 0.5) Explicitly, ao reaches 0.5 at a s ss 2 using LO evolution for ao, but a s « 1.3 using NLO. It is encouraging 
that the trend to smaller (more negative) AG at low Q 2 is enhanced (rather than cancelled, as it certainly might have 
been) by considering the NLO corrections of Eq. (||) . 

Moving the other direction, one can ask how far the trend of increasing a s AG with increasing Q 2 continues. From 
Eq. ([jj), as Q 2 increases, ao(Q 2 ) (slowly) decreases and this in turn reduces the rate of change of a s AG. In this way, 
a s AG inevitably increases, but assuming that AE is in the range (.5 <-> 1), Eq. ( |l0| ) shows that if a s AG reaches 
around (1 <-> 2), the increase with scale effectively halts, and a AG becomes renormalization group invariant beyond 
second order in a s , as well. 

Following the logic of refs (9j and JTTj[ ] , we might argue that evolution of both unpolarized momentum fractions and 
polarized spin fractions are consistent with a nonrclativistic quark picture at a low energy scale of a s (/i§)/7r w 1±0.2, 
where the gluonic contributions vanish. Thus, with AG(/i§) = 0, and AE = .5 p0| , we would predict, to leading 
non-zero order, 

a s AG(5 GeV 2 ) = - (4//3 )(a s (5 GeV 2 ) - a s (/^ 2 ))(AE) « 0.6 ± .1, (12) 

i.e. AG(5 GeV 2 ) « 2.1 ± .4, a result which is certainly compatible with recent high energy experimental estimates 
, even though we begin with AG = at the low energy scale. 

Jaffe provides us with more sophisticated model predictions for AG at the low energy scale ||, namely for the 
non-relativistic quark model AG(NRQM) = -0.8 and a s (NRQM) = 0.9, while for the bag model AG(bag) = -0.2, 
and a s (bag) = 2.0. His paper makes it clear that these numbers are not to be taken too seriously, but the sign of AG is 
taken as an interesting and surprising result. However, if we evolve these as above (bearing in mind the caveats outlined 
in the introduction) we would predict a s AG(5 GeV 2 )(from NRQM) = -.5, while a s AG(5 GeV 2 )(from bag) = +.1. 
The latter has shown a sign change, with AG(5 GeV 2 ) = +.35 only slightly below the low end of the experimental 
range. Assuming a (scale independent) value of AE = .7, rather than .5, brings these predictions even higher, e.g. 
from the bag model result it predicts AG(5 GeV 2 ) = 1. We see that the sign can indeed change, although it depends 
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on both the starting values for a s and AG, and thus on the quantitative details of the model. This is also clear from 
the evolution shown in Figs. 2 and 3, where the location (in a s ) of the sign change is a strong function of the starting 
values. 

We conclude that the sign of a s AG at low energy scales is ill-determined simply based on 0(a 2 ) (which is leading 
non- vanishing order) perturbative QCD arguments, starting from experimental values at moderate Q 2 . Both LO and 
NLO perturbative evolution tend to decrease a s AG, even to negative values, depending on starting conditions. The 
spin content of the glue remains an important aspect of nucleon structure, one whose measurements should certainly 
be improved at high Q 2 , and which may play a key role in understanding and interpreting low energy models, but the 
connection between the two scales involves nonperturbative physics and is thus not trivial to predict, and one should 
not be particularly surprised if aAG is found to be negative in a quark model calculation. 
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FIG. 1. Plot of the singlet first moment ao(Q 2 ), as a function of a s /iv. Solid (dashed) curve shows leading (next to leading) 
order calculations. We started from experimental values Q at 5 GeV 2 and evolve downwards in Q 2 . The two lines represent a 
band of possible results, beginning at ao(Q 2 = 5 GeV 2 ) = .37 ± .11. The horizontal dotted line shows ao = 0.6, which would 
be the rough expectation in the naive quark model picture, assuming that AE (scale independent) corresponds to the "spin on 
the quarks" , when there is no significant contribution from the gluon anomaly term in Eq. (^) . 
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FIG. 2. Same as previous figure, but showing a s AG, as calculated in Eq. (|lo|), as a function of a s /w. The dashed curve 
corresponds to using ao as computed in Eq. (^), to NLO. The horizontal dotted line is a s AG = 0, where the gluon anomaly 
term makes no contribution to Eq. m). 
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FIG. 3. Same as previous figure, but showing just AG as a function of a 3 /w, evaluated using Eq. ([lo|). The dashed 
corresponds to using ag as computed in Eq. fel to NLO. 



7 



